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ABSTRACT 

Spectroscopic orbits have been reported for nine unseen companions orbiting solar- 
type stars with minimum possible masses in the range 0.5 to 10 Jupiter masses. We 
compare the mass distribution of these nine planet candidates with the distribution 
of low-mass secondaries in spectroscopic binaries. Although we still have only a very 
small number of systems, the two distributions suggest two distinctive populations. The 
transition region between the two populations might be at the range of 10-30 Jupiter 
masses. 

Subject headings: binaries: spectroscopic — planetary systems 
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1. INTRODUCTION 

Eight candidates for extrasolar planets have been 
announced over the past two years (e.g.. Marcy & 
Butler 1998). In each case, very precise stellar radial- 
velocity measurements, with a precision of about 
10 m s^^ or better, indicated the presence of a low- 
mass unseen companion orbiting a nearby solar-type 
star. The individual masses of the eight companions 
are not known, because the inclination angles of their 
orbital planes relative to our line of sight could not 
have been measured. The minimum masses for the 
eight candidates, attained for an inclination angle of 
90°, are in the range 0.5 to 7.4 Jupiter masses (Mj). 
These findings render the eight companions to be gi- 
ant planets or at least 'planet candidates'. 

The detections of these eight companions were an- 
nounced seven to nine years after a companion of HD 
114762 was discovered (Latham et al. 1989), based on 
measurements with a lower precision (Latham 1985). 
Mazeh, Latham, & Stefanik (1996) have shown that 
the minimum mass for the companion of HD 114762 
is 9.4 Mj. Therefore, when considering the emerging 
population of planet candidates, HD 114762 should 
be considered together with the eight new candidates. 
Table 1 lists the minimum mass, period and discovery 
date of the nine objects. For random orbital orienta- 
tions, the expectation value for sin i is 0.76, so the 
actual masses of the nine companions are expected to 
be in the range of 0.6-12 Mj. 

The nature of the newly discovered low-mass com- 
panions is not yet clear. They could be planets, as 
suggested by various authors (e.g., Marcy & But- 
ler 1998), or many could just be brown-dwarf secon- 
daries, formed in binary stars (Black 1997). With the 
small, but not insignificant number of spectroscopic 
orbits implying planetary minimum masses, we can 
now begin to study the distribution of their orbital 
parameters, in order to address this very basic ques- 
tion. 

In this paper we discuss the emerging difference be- 
tween the mass distribution of planet candidates and 
the low-mass end of the distribution of binary secon- 
daries. This point has been already discussed by pre- 
vious studies (Basri & Marcy 1997; Mayor, Queloz 
& Udry 1998; Mayor, Udry & Queloz 1998; Marcy 
& Butler 1998), but in those papers the mass distri- 
bution was binned linearly. Here we choose to use 
a logarithmic scale to study the mass distribution, 
because of the large range of masses, 0.5-300 Mj, in- 



volved. The logarithmic scale has also been used by 
Tokovinin (1992) to study the secondary mass distri- 
bution in spectroscopic binaries, and was suggested 
by Black (1998) to study the mass distribution of the 
planetary-mass companions. 

This work is based on an extremely small sample, 
and the validity of our results will need to be verified 
by many more detections. However, if verified, the 
difference in mass distributions that we find might 
provide an important clue for how to distinguish be- 
tween planets and low-mass stellar companions. 

A preliminary version of this work was presented 
at the meeting "Physical Processes in Astrophysical 
Fluids", in Haifa, January 1998 (Mazeh 1998). 

2. COMBINED MASS DISTRIBUTION 

We wish to consider the mass distribution of the 
planet candidates and compare it with that of the 
low-mass secondaries in spectroscopic binaries. To do 
so we use the results of two very large radial-velocity 
studies of spectroscopic binaries recx^ntly completed, 
for which partial results have been published. One 
sample is composed of G and K stars studied by 
Mayor et al. (1997), and the other is the Car- 
ney & Latham (1987, hereafter C-L) high-proper- 
motion sample. We will use the first sample to 
estimate the mass distribution in the mass range 

1 < log{M/Mj) < 2, and the latter sample to es- 
timate the value of the distribution in the range 

2 < log{M/Mj) < 2.5. 

Mayor et al. (1997) listed 10 spectroscopic binaries 
with minimum secondary masses in the range 10-63 
Mj. Their list includes only two systems in the range 
1 < log{M/Mj) < 1.5 and 8 systems in the range 
1.5 < log{M/Mj) < 2. We note that their list does 
not cover the second range completely, because the 
table of Mayor et al. does not include binaries with 
minimum secondary masses between 63 and 100 Mj. 
Mayor et al. were kind enough to let us know that 
they have found 5 additional binaries in the range of 
63-100 Mj (Halbwachs, private communication). 

Detailed results for the C-L sample are not yet pub- 
lished, but partial results were presented in two con- 
ference papers (Latham et al. 1998; Mazeh, Goldberg 
and Latham 1998). Mazeh, Goldberg and Latham di- 
vided the sample into two subsamples, with high- and 
low-mass primaries. We use here only the high-mass 
primary subsample, with primary masses between 0.7 
and 0.85 Mq, because they are more similar to the 
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Table 1 
The Planet-Candidates 



Name 


M2 sin i 


P 


Discovery Ref. 




(Mj) 


(days) 


Date 



HD 114762 


9.4 


84 


1989 


1,2 


51 Peg 


0.5 


4.2 


1995 


3 


47 UMa 


2.5 


1090 


1996 


4 


70 Vir 


7.4 


117 


1996 


5 


55 Cnc 


0.8 


14.7 


1996 


6 


T Boo 


3.9 


3.3 


1996 


6 


V And 


0.7 


4.6 


1996 


6 


16 Cyg B 


1.6 


804 


1996 


7 


p CrB 


1.1 


39.6 


1997 


8 



Note. — ^Latham ct al. 1989; ^Mazeli, Latham 
& Stcfanik 1996; ^Mayor & Qucloz 1995; '^Butler & 
Marcy 1996; ^Marcy & Butler 1996; ^Butler et al. 
1997; ^Cochran et al. 1997; ^Noyes et al. 1997. 



primaries in the other samples considered here. From 
Figure 1 of Mazeh, Goldberg and Latham we can es- 
timate the number of systems with secondary masses 
in the range 2 < log{M/Mj) < 2.5 to be 20. 

The results of the two samples of spectroscopic bi- 
naries have to be scaled to the size of the sample out 
of which the nine planet candidates were found. The 
scaling is not simple because the parent samples in 
which the planets were found are not well defined. 
The nine planets were discovered by different research 
groups, with different time coverage and slightly dif- 
ferent precision (e.g., Marcy & Butler 1998). For the 
present discussion we will assume that the total num- 
ber of stars searched was two hundred. This has to be 
compared to the 570 stars of Mayor et al. (1997) and 
420 stars of the sub-sample of high-mass primaries of 
the C-L sample (Mazeh, Goldberg & Latham 1998). 
The results are summarized in Table 2, where Nsd 
is the number of binaries, scaled to a sample of 200 
systems. The combined scaled histogram is plotted 
in Figure 1. 



0.5 1 1.5 

Secondary Mass [log (M /M. )] 



Fig. 1. Scaled histogram of the cxtrasolar planet- 
candidates and the low-mass secondaries of spec- 
troscopic binaries. The dashed line is the stel- 
lar /substellar limit. 
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Table 2 

The Scaled Mass Distribution 



Mass Range 


# of Observed Systems 


Scaling 


Nscl 


(Mj) 




Factor 




-0.5 < log{M) < 0. 


3 




3 


0. < log{M) < 0.5 


3 




3 


0.5 < logiM) < 1. 


3 




3 


1. < log{M) < 1.5 


2 


200/570 


0.7 


1.5 < log{M) < 2. 


8 + 5t 


200/570 


4.6 


2. < log{M) < 2.5 


20 


200/420 


9.5 



Note. — t unpublished data (Halbwachs, private communication) 



3. ESTIMATE OF THE CORRECTED DIS- 
TRIBUTION 

Before considering the possible interpretation of 
the combined histogram, we have to correct the his- 
togram for two effects. The first one has to do with 
the fact that the masses given in Table 1 and in Mayor 
et al. (1997) list are only minimum masses, and there- 
fore the actual mass of each secondary is most proba- 
bly larger. The correction of this cflFect tends to shift 
the distribution towards the right side of Figure 1. 
The second effect reflects the fact that binaries with 
too small amplitudes could not have been detected, 
because their period is too large, or their inclination 
angle is too small. The correction of this effect tends 
to increase the number of companions detected in bins 
with small masses, while the effect is negligible for 
bins with large masses. Both effects were taken into 
account in the work of Mazeh, Goldberg & Latham 
(1998), so we need to correct only the counts of the 
two other samples. 

To correct for the first effect we calculated the 
probability of each system to fall in each bin of the 
histogram, assuming random orientation in space. To 
derive a modified histogram we added up the contri- 
butions of each binary to each bin of the histogram, 
denoting the resulting counts by A^mod- 

To correct for the second effect we consider the 
probability of not detecting a binary or a planet 



in a systematic radial-velocity search (see Mazeh, 
Latham & Stefanik (1996) for details). Suppose that 
the search detects all stars with radial-velocity mod- 
ulation with a period P between Pmin and Pmax, 
and with semi-amplitude K larger than or equal to 
the search threshold Kj^in- For given primary and 
secondary masses, and for a given orbital period, 
all systems with an inclination smaller than some 
threshold inclination cannot be detected, because K 
is smaller than JsTmin- We can therefore calculate 
U[K^in]{P, Ml, M2) — the probability of not detect- 
ing a binary by a search with a given threshold Kj^in, 
assuming random orientation in space (e.g., Mazeh & 
Goldberg 1992). 

To get the probability of not detecting a binary 
taken at random from a population of binaries, with 
given range of secondary masses and periods, we 
have to integrate U[Kjnin]{P, Mi, M2) over the given 
parameters of the population. We will then get 
U[K„i[n]{Mi), which presents the probability of not 
detecting a binary, averaged over the secondary masses 
and period domains. 

To correct for the undetected binaries we have to 
multiply the number of systems in each bin by the 
corresponding 

c=(i-u[K^i4y\ (1) 

where we have dropped the dependence on Mi. 
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The main parameter here is Kmin- wliic;li in tm^n 
strongly depends on the precision per measurement, 
but also on the number of measurements per star and 
their temporal distribution. Therefore, the exact val- 
ues of i^Tmin for each of the samples discussed here are 
still not well known. For the planet searches we will 
assume that Xmin is twice the precision per measure- 
ment. This means that any binary with a peak-to- 
peak variation larger than 2Kjnin, or four times larger 
than the error per measurement of the survey, was de- 
tected. We therefore assume iCmin to be 20 m s~^ for 
the planet search. For the Mayor et al. sample the 
detection threshold is larger than twice the precision 
per measurements (Halbwachs, private communica- 
tion), so we will assume, somewhat arbitrarily, Kmin 
of 1 km s~^. 

To calculate the correction factor for each bin of 
the histogram, we considered a population of binaries 
with secondary mass range coinciding with the bin 
mass range, with a Duquennoy & Mayor (1991) pe- 
riod distribution between 1 and 1500 days. We then 
applied the correction derived to the modified counts 
in each bin to get Ncor- We estimated the error in 
each of the first three bins by the square root of the 
modified number of systems in each bin, multiplied 
by the correction factor. For the Mayor et al. sample 
we took into account the fact that the original sample 
was larger by 570/200, so the relative errors for these 
two bins are smaller by the square root of this scaling 
factor. The 'corrected' histogram is given in Table 3. 

Note that both samples cover the range 1. < 
log(M / M j) < 1.5, both of which yielded very similar 
estimates. When plotting the corrected histogram in 
Figure 2 we combined the two estimates together and 
got 1.5 ± 1.0 for this bin. 

4. DISCUSSION 

The corrected combined histogram might suggest 
that we sec here two populations. At the high-mass 
end of the histogram we see a distribution which drops 
steeply when we move from 200 to 20 Mj. At the 
planetary range of masses we see a flat distribution, 
which might even rise very mildly when we move from, 
say, 20 to 0.6 Mj. Unfortunately, the number of sys- 
tems in each bin is small. However, the two different 
slopes in the two parts of the diagram seem real, as 
they are based on more than one bin. 

The derived diagram depends on two parameters 
— Kmin and the number of bins of the histogram. 



We got the same gross featiires, namely two oppo- 
site slopes in the two parts of the diagram, when we 
changed the values of these two parameters. Dividing 
the total range of the diagram into 5 or 4 bins instead 
of 6 bins shifted the transition region between the two 
slopes somewhat to the right. Changing K^in from 
20 to 50 m s^^ made the slope at the left hand side of 
the diagram steeper. We conclude therefore that the 
overall shape of the diagram does not depend strongly 
on the specific values of the parameters of the deriva- 
tion. 

The transition region between the two populations 
is at about 10-30 Mj. Unfortunately, the relative er- 
ror of this bin is very large. Nevertheless, it seems 
that this is the bin with the smallest number of sys- 
tems. The very low comit estimate in this bin is sup- 
ported by the fact that the very sensitive searches for 
planets, which yielded the discovery of the eight new 
planet candidates, did not find any companions with 
minimum masses between 10 and 30 Mj. With i^min 
of 20 m s~^ these searches could detect more than 
99% of the binaries in this bin. The fact that no pre- 
cise search discovered any binary in this bin indicates 
that the number of systems with secondary masses 
between 10 and 30 Mj is very small. 

The drop of the secondary mass distribution we 
find when moving from 200 to 20 Mj is consistent 
with the finding of Halbwachs, Mayor & Udry (1998), 
who studied the mass ratio distribution of spectro- 
scopic binaries in the samples of G and K stars of 
Mayor et al. (1997). Halbwachs, Mayor & Udry 
found a flat histogram of the mass ratio, although 
they could not exclude increasing or decreasing power 
laws of the form q", where q is the mass ratio and 
-0.82 < a < 0.87. The flat distribution of q yields 
constant dN/dm^, if all primary masses are similar. 
This corresponds to dN/d\og{m2) oc m2, consistent 
with our findings. The drop we find is also consistent 
with the findings of Mayor, Queloz & Udry (1998; 
see also Mayor, Udry & Queloz 1998) who found 
that dN/dm2 cx: m^*^'^. Their result corresponds to 
dN/d\og{m2) oc mj°'^. Figure 2 of this work sug- 
gests a steeper drop, but the difference is within the 
errors. 

The transition region between the two populations, 
or between the two slopes, that we find here is, how- 
ever, different from the findings of Mayor, Queloz 
& Udry (1998) and Mayor, Udry & Queloz (1998). 
They find a borderline at 7 Mj, while our logarith- 
mic treatment of the data suggests a transition region 
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Table 3 

Estimated Corrected Mass Distribution 



Mass Range Nmod Correction 
(Mj) Factor 


N 

cor 


Planet Search 


-0.5 < log{M) < 0. 1.6 2.2 
0. < logiM) < 0.5 2.7 1.1 

0. 5 < log{M) < 1. 2.6 1 

1. < logiM) < 1.5 1.5 1 


3.6 ±2.8 
2.9 ± 1.8 
2.6 ± 1.6 
1.5 ± 1.2 


Mayor ct al. Sample 


1. < log{M) < 1.5 0.4 3.9 
1.5 < log{M) < 2. 3.1 1.3 


1.6 ± 1.5 
4.0 ± 1.4 


Carney & Latham Sample 


2. < log{M) < 2.5 


9.5 ±2.1 



10 - 



8 - 
6 - 
4 - 



2 - 



-0.5 0.5 1 1.5 2 2.5 

Secondary Mass [log (Mg^^/M.^p)] 

Fig. 2. — Corrected histogram of the extrasolar 
planet-candidates and the low-mass secondaries of 
spectroscopic binaries. The dashed line is the stel- 
lar /substellar limit. 



at the range of 10-30 Mj. Another difference is the 

shape of the distribution in the planetary mass range. 
They find a very steep rising distribution when mov- 
ing down towards the range of 1-5 Mj. We find an 
almost flat logarithmic distribution, with perhaps a 
mild rise towards lower masses, depending on the ex- 
act value of /Crnin- 

Let us assume that Figure 2 indeed shows two dis- 
tinctive slopes in the two parts of the diagram. Let 
us further assume that this reflects the fact that we 
sec here two different populations, one below 10 30 
Mj, and one with masses larger than this transition 
region. One possible interpretation of the diagram, if 
indeed we see here two different populations, is that 
the two populations were formed differently. Maybe 
the lower-mass population was formed like planets, 
out of an accretion disc, while the higher-mass pop- 
ulation was formed like binary stars, in a mechanism 
which probably involves large-scale gravitational col- 
lapse (e.g. Boss 1996; Black 1986). If this is the case 
then the binary secondaries include stars and brown 
dwarfs together. 

Figure 2 suggests that the transition region be- 
tween the two populations is at about 10-30 Mj. 
This is of astrophysical significance, if indeed the 
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lowor-mass population is composed of planets, as it 
might tell us about the lower limit and upper limit 
of the formation of secondaries and planets, respec- 
tively (Marcy & Butler 1995, 1998; Mayor, Queloz & 
Udry 1998; Mayor, Udry & Queloz 1998). The upper 
limit of the planetary masses is set by the conditions 
in the accretion disc, and most probably by the in- 
teraction between the planet and the gas and dust in 
the disc. Boss (1996) already noted that Lin and Pa- 
paloizou (1980) theoretically predicted that the max- 
imum mass for the formation of a planet in the disc 
of the Solar nebula is about 1 Mj. As the maxi- 
mum mass depends on the mass of the early neb- 
ula, we can get somewhat higher masses in different 
cases. The lower limit for secondary masses in bina- 
ries is set by the binary formation mechanism, what- 
ever that mechanism might be. Boss (1988), for ex- 
ample, noted that the theory of opacity limited cloud 
fragmentation predicts that the minimum mass for a 
companion is about 10 Mj. In fact. Low & Lynden- 
Bell (1976) estimated already 20 years ago that the 
minimum Jeans mass for fragmentation of a molecu- 
lar cloud is 7 Mj. Silk (1977), in a contemporaneous 
study, came up with minimum masses between 10 and 
100 Mj, depending on the shape of the collapse. If 
we indeed see the transition region between the two 
populations at about 10-30 Mj, this is not too far 
from the predictions of the theories. 

Duquennoy & Mayor (1991; see also Mayor, Queloz 
& Udry 1998) have suggested that the observed or- 
bital eccentricities can be used to distinguish between 
planets and stellar companions. However, Mazeh, 
Mayor & Latham (1996), when discussing the eccen- 
tricity versus mass of the known planet candidates, 
pointed out that the planet-disc interaction (e.g., Gol- 
dreich & Tremaine 1980) is a possible mechanism for 
generating a strong dependence of eccentricity versus 
mass (Artymowicz 1992; Lubow & Artymowicz 1996), 
at least for moderate eccentricities. This possibility 
can undermine the potential of the eccentricity-mass 
dependence to distinguish between planets and sec- 
ondaries. Furthermore, Black (1997) analyzed the ec- 
centricity as a function of period and concluded that 
the eccentricities observed are consistent with the as- 
sumption that all the planet candidates are actually 
low-mass brown dwarfs formed like binary stars. It 
seems therefore that it might be too early to distin- 
guish between brown dwarfs and planets solely on the 
basis of their orbital eccentricity. 

Mazeh, Mayor & Latham (1996) speculated that 



"The 10 40 Mj mass gap may prove to be critical 
for the interpretation of" the eccentricity-mass de- 
pendence. We confirm here that the transition region 
between the two populations could be at this range 
of masses. 

Obviously, the left hand side of the histogram and 
the transition region between the two slopes derived 
in this paper are based on a very small number of 
objects all together, and these features need to be 
verified by many more detections. Further, one still 
needs to make sure that the different slope in the 
planetary-mass range is not due to some selection 
effects. For example, there might be a correlation 
between the orbital period and the secondary mass, 
which might make the small-mass secondaries easier 
to detect. Such an effect could cause the histogram 
to appear to rise towards smaller mass. However, if 
the shape of the histogram can be verified, and if the 
planetary-mass objects prove to be extrasolar plan- 
ets, the shape of the histogram might give us the 
long-sought clue for how to distinguish planets from 
low-mass stellar companions. 

We express our thanks to J.-L. Halbwachs for his 
very useful comments on the manuscript. We thank 
M. Mayor, S. Udry and J.-L. Halbwachs for letting us 
use their unpublished results. We thank the referee, 
Dr. D. Black, for a critical reading of the manuscript 
and for his comments that led to significant improve- 
ment of the paper. This work was supported by US- 
Israel Binational Science Foundation grant 94-00284 
and by the Israeli Science Foundation. 
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